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Abstract 

The D — > X-ul^l^ transitions - branching ratios, forward-backward asymmetry Apg, the CP asymmetry 
and the CP asymmetry in the forward-backward asymmetry - have two sources: for D± they 
represent a pure AC = 1 & AQ = current interaction whereas neutral D mesons can also communicate 
via their anti-hadron. Standard Model (SM) contributions to BR{D — >■ Xul~^l~) come primarily from long 
distance dynamics, which overshadow short distance contributions by several orders of magnitude; still they 
fall much below the present upper experimental bounds. Even the SM contributions to Ap^, A^p and Apg 
are tiny, quite unlike in beauty hadrons. The branching ratios are hardly dented by contributions from 
the Littlest Higgs Models with T parity (LHT) even in the short distance regime, let alone in the SM long 
distances dynamics. Yet the asymmetries Ap^, A^p and Ap^ in these New Physics models can be enhanced 
over SM predictions, as they arise purely from short distance dynamics; this can occur in particular for 
ApQ and Ap^ which get enhanced by orders of magnitudes. Even such enhancements hardly reach absolute 
sizes for observable experimental effects for Ap^ and A'^p. However LHT contributions to Ap^ could be 
measured in experiments like the LHCb and the SuperB Collaboration. These results lead us to draw further 
conclusions on FCNCs within LHT-like models through some simple scaling arguments that encapsulate the 
essence of flavour dynamics in and beyond the Standard Model. 
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1 Introduction 

The discovery of charm quarks was and still is seen as a great success of the Standard Model (SM), since 
their existence was necessary for the observed suppression of strangeness changing neutral currents. Charm 
hadrons were also found in the expected mass range as was the predicted preference for decays to strange 
hadrons. 

Yet, at last, some possible hint of New Physics (NP) has appeared in charm physics. Compelling evidence 
for — Z)° oscillations has been presented by Belle, BaBar and CDF [1]. The HFAG has combined the 
results on neutral D decays allowing CP violation [2]^: 

.z. = ^ 0.63t°- , ,z, = 1^ = 0.75 ±0.12 

0.9lt°:l« , = -lO.2t^:^n (1) 

The observation of — oscillations is hardly disputed, while the relative size of xd and t/^i is not clear 
yet. Before these experimental results, most theorists argued that the SM predicts xd, yo < 10""* - yet 
not all: in 1998, xd, Vd < 10~^ was called a SM conservative bound [3]; in 2000 and 2003 a SM prediction 



^Up to date results can be found in the HFAG website 
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2 SM Contributions to D ^ XJ+T 



obtained from a sophisticated operator product expansion yielded xd, Ud ^ 0(10^'^) [4] and more recently 
in [5]; alternatively, in 2001 and 2004 a SM prediction on D° — I)° oscillations was based on SU (3) breaking 
mostly in the phase space for yu and then from a dispersion relation for xd [6]. 

While the present experimental results on xd and yo can be accommodated within the available theoret- 
ical SM estimates, and no non-zero CP asymmetry has been seen yet, the observation of — oscillation, 
however, has 'wetted' the appetite on thinking of NP in charm decays. The authors of Ref. [7] consider 
a (approximately) S'?7(2)i-invariant NP scenario; therefore NP contributions to — and — 
oscillations are not independent of each other. 

There is a large variety of NP models, in which the Higgs boson appears as a pseudo-Nambu-Goldstone 
boson of a spontaneously broken global symmetry, namely the 'Little Higgs' class of models [8,9]^. To 
achieve this program one needs at least heavy gauge bosons Zn and A//, a heavy top partner T and 
a scalar triplet $ as physical degrees of freedom as is implemented in the 'Littlest Higgs' model [11,12]. 
Studying electroweak precision observables shows that for such new states to arise below the 1 TeV scale, 
one needs an additional discrete symmetry [13, 14], called T parity: the SM particles and the heavy top 
partner T are even and M^^, Zh^ -Ah and $ are odd. A consistent implementation of T parity requires also 
the introduction of the so-called 'mirror' fcrmions - one for each quark and lepton species - that are odd 
under T parity [14,15]. This creates the 'Littlest Higgs' model with T parity (LHT)'^. While some theorists 
probably see it as intellectually economical, most experimentalists do not view it like that; at the same 
time they should understand that it can provide them with more work, but less so than SUSY! The most 
important point here is that the motivation for LHT models come from outside flavour dynamics; at the same 
time they can create important non-trivial signals of NP in B, K - and D physics. The LHT models can - 
not necessarily, but possibly - affect AC — 2 dynamics significantly. In particularly, it can generate sizable 
or even relatively large indirect CP violation in D*^ decays [17]. It can implement a dynamical realization of 
the symmetry approach described in Ref. [7]. 

Encouraged by the findings of Ref. [17] about the possible impacts of LHT models in AC ~ 2 dynamics 
we had looked at two AC = 1 processes, Z?" 77 and D° — >■ in a previous study [18]. While LHT 

failed to contribute significantly to the total decay rates in these channels, it had led us to some general 
conjectures on charm changing neutral currents (CCNC) within LHT-like scenarios. In this current work we 
continue to analyze the impact of LHT on another AC = 1 process to probe deeper into CCNCs in LHT-like 
scenarios, namely to D — >■ Xul~^l~. Similar to our previous work, we do not see any sizable enhancements 
in the global decay rates. However in the presence of large weak phases in LHT-like models it might seem 
surprising to find also very tiny contributions to the CP and forward-backward asymmetries and . 
This is because the (non-trivial) SM asymmetries get produced by short distance dynamics, and LHT-like 
scenarios can create much larger A^p and Apg than the SM can; yet they are still small in their absolute 
size. However the LHT contributions to the forward-backward CP asymmetry Apg are sizable - even large 
- such that they can be experimentally measured in the coming decade. 

In this article, we discuss both short and long distance SM contributions to D ^ Xul^l^ in Sect. 2 along 
with Aqp, ApQ and ^pg in this channel. We will describe mostly Xul^l^, since it is given only 

by AC = 1 couplings, while one can also produce D° Xul~^l~ by AC = 2 couplings due to — 
oscillations. We go on to briefly introduce LHT and its contributions to Z) — > X^l^l^ in Sect. 3. Our 
quantitative findings are presented in Sect. 4. We take a more critical look at CCNCs in LHT-like models in 
Sect. 5. In Sect. 6 we put forward some simple scaling arguments to explain why NP interventions such as in 
LHT-like models have effects of the size we see and the conclusion to this work follows in Sect. 7. 

2 SM Contributions to D ^ XJ+l' 

The transition of Z? — !■ Xul^l^ must be produced by charm changing neutral currents, which are much weaker 
even than their strangeness and beauty analogues in the SM. These decay rates are tiny and dominated by 
long distances effects. Yet the forward-backward Apg, CP asymmetries A^p and CP asymmetry in the 
forward-backward asymmetry Ap^ could still be controlled by SM short distances dynamics. However 
the dynamics scenery is very complex as shown below; conceptually very similar to _B — > Xsl^l~ , but 

^For an overview of the different 'flavours' of Little Higgs models cf. [10] 
^For a detailed description of tlio Littlest Higgs Model with T parity cf. [16]. 
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2 SM Contributions to D ^ XJ+T 



quantitatively at a smaller level. One can learn a lot by studying these asymmetries about SD dynamics 
in general where different operators mix and an alternative perspective can be obtained from B studies. It 
turns out - not surprisingly - that a careful scrutiny needs huge statistics. This would be an important task 
for a Supcr-B factory like the recently approved SuperB project undertaken by the INFN [19]; LHCb might 
also be able to address it. 

We will discuss first Xul^l^, since it proceeds purely by a AC = 1 interaction; then we will 

comments on lessons learnt from neutral D — > transitions, where — oscillations can get 

involved. 

2.1 Tsm{D XJ+1-) 

The quark level process c — >■ ul^l^ is described with an operator basis of the following ten operators: 



05 = (uM)^(g-^,7^<?^,) 



1 

06 = «7mC^)E(9^.T^'?S,) 



O7 
O9 



167r2 



■{uLj^,CL)ih''l) Oio = 



167r2 

^2 



T{uLlt,CL)(ll^lhl) 



(2) 



where q = d, s,b and a, /3 are colour indices. The charm operators Oi, i ~ 1, 10 are analogous to those in 
b decays [21]. The effective weak Hamiltonian is expressed in terms of these operators taken at a scale fi: 



i=3 



(3) 



where the coefficients describe renormalization of the operators at /x from the normalization scale. For the 
starting point one can naturally chose Mw'- 

c[''\Mw)^0 C^«)(Mw) = 1 C3-6(A/w) = 

j=d,s 



j=d,s 



j=d,s 



j=d,s 



sin {9w) 
Yojxj) ~ Yo{xb) 



- 4Co(xj) - Do{xj) - {xj Xb) 



(4) 



where Xj = m'j/m'^, with rrij being the masses of the internal down type quarks j = d,s,b. Here we 
have used the unitarity of the CKM matrix, J2j=d s b ^uj ^jc = to eliminate the dependence of the Wilson 
coefficients on the third family. The limit — > cannot be taken as there is a logarithmic divergence in 
D^{x) in that limit. The form factors arc defined in a modified way from Ref. [22]: 
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2 SM Contributions to D ^ XJ+T 



Coix) 



1 /x 3 1 

2 \ 4 ^ 



3 2x2 _ ^ 



Do{x) = QdEo{x 



(9 - 43x + 28^2) X (21 - 66a; + 41a;2 _ 2x^) 



12(1 -a;)3 



2 , , x2 (15- 16a; + 4x2) 
3^°s(-) + 6(T3^)^ ^°S(- 



12(1 - x)4 
X (18 - llx - x^l 



log(x) - 27(x) 



D'oix) = QdE'^{x) 



12(1 -x)3 

1 11x2 _ 7^, + 2 6x3 

3+ 4(1 -x)3 +4(r^^°g^") 



£;^(x) 



5 

12 



1 - 5x - 2x2 



3x2 



4(x-l)3 ' 2(x-l)4^°S(^) 



3x 



X — 1 (x — 1)' 



■ log(x) 



(5) 



Qd is the charge of the internal down type quarks. Here 7(x) is the gauge dependent term which, for = 1, 

4 



IS 



7(x) 



8 V(a^-l)' 



■ log(x) 



X- 1 



(6) 



At scales /i < Mw one can express HeffifJ.) using Ci{fi), i = 1,...,10 evolving through the two-foop QCD 
renormaUzation group equation. Actually one has two regimes, Mw to rrib and rrif, to mj'] a matching con- 
dition as usual is applied through as(mb, mb; 4) = ^.^(mb, mb; 5). These operators mix via rcnormalization; 
in particular O7 mixes with Oi, O2 [21,23], O^-q and 0%. However, as pointed out in [23], O7 is completely 
dominated by the two loop QCD radiative correction, which was taken into account in [24,25]: 



as{nic) 



8 6 



i=l j=l 



An 



-C2(mc) 



where 



Vb = -, T' ^ 



(7) 



o.s{mc) 



(8) 



The matrix X, the vector z and the function / introduced above are given in Appendix A along with the 
Wilson coefficients for the operators Oi_6 both at /i = rrih and /i ~ rric- We include the next to leading order 
correction to the running of the as with as{mw) = 0.125, as(TOt) = Q!s(TOh, mw , 5), ag{mc) = as{mc^ mb, 4). 



v{fi,fj.',nf) 
/3i('V) 



as{n',fi',nf) / _ f3i{nf) Qs(/i^^^n/) log(t>(^, ^^ n/)) 
v{fi,fi',nf) V Po{nf) 47r v{fi,n',nf) 

I- /^o(.v)^^^^4^1og(^) 

2 

II - 3^/ 
102 - 



(9) 



Although some of the formfactors have gauge dependence, it drops out of the final amplitude as expected. 
^Of course, one can question the robustness of the selection of ^ = mj, and fj, = rric- However, any arguments on this choice 
arc purely acedemic at this point considering the level of precision aimed at in the current calculations. 
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2 SM Contributions to D ^ XJ+T 



Og mixes with Oi-e beyond the leading order. However, only Ci and C2 are numerically significant and 
lead to an important cancellation amongst themselves. As we shall later see, the dominant contribution to 
the SM branching fraction comes from Og; hence it is important to have a detailed look at it contrary to 
what was argued in [24,25]. The Wilson coefficient for Og after QCD corrections is given by ^ 



Cg(mc) = Csimw) + K^cjhiz^^s) ^ 3a(m,) + C,+i(7«,) (10) 

j=d,s y=l,3,5 J 

where h{zj, s) comes from the one loop QCD correction to the four fermion operator and is given by ^ 
h{z,s) ^ Qdh{z,s) 



7, 4, mc 4, , , 4 8z2 1 / 4^2^ 

h{z,s) = --log log 2 +- + ---- 2 + — 

3 fi 3 93s3\ s 



1 



4z2 



2 tan 



-1 1 



log 



if s < 4z^ 

in if s > 4z^ 
(11) 



with 

. _ (P;+ + Pi- )^ _ 
s — „ , — 

rric 

As in the case of the analogous decay in B systems, the logarithmic term in h{z, s) exactly cancels the 
logarithmic dependence in Cg{m\Y) that comes from the electromagnetic penguin and hence removes the 
logarithmic dependence on light quark masses at the scale n = mw as was pointed out in [24]. The 
importance of QCD correction was pointed out in [27], but we disagree with their argument that the purely 
electroweak "Inami-Lim" contribution to Cg should be ignored as it is dependent on light quark mass and 
is reproduced as a limit of the QCD correction when s — > 0. A careful look at the form of /i(z, s) shows that 
the logarithms have opposite signs in the "Inami-Lim" term and in the QCD correction as argued before*. A 
discussion of the logarithmic dependence of Cg before including QCD corrections and its cancellation after 
including the same is discussed in [28] for the case of B mesons. A similar argument applies in the case of 
the D mesons too. 

Oio does not suffer from any QCD corrections [24, 26] which makes the assumption made in [25] unnec- 
essary ^ . Finally, the differential decay branching fraction is given by 



-^Br^^ {D ^ XJ+l 
1 Glo 



D 



7687r5 



■-{l-sY 



{\C^{^l)f + \CMf) (1 + 2s) + 12 Re(C7(/x)Cg*(/i)) + 4 



1 + ?)|C7(m)|' 



(12) 



with fi = rric = 1.2 GeV. Integrating over s gives us the total decay rate. One has to be careful about not 
picking up the infrared divergence in the differential decay rate. We made an infrared cut on s at about an 
invariant dilepton momentum of 20 MeV. We get a branching fraction of 



BRl^{D XuC+e-) - 3.7 x 10"^ 



(13) 



which is smaller than what is stated in [24, 25, 29] but larger than the number in [27] for reasons stated 
above. The same for muons in the final state is slightly smaller due to the finite mass of the muon. The SD 



^The dependence of C<){mc) on h{l,s) and h{0,s) have been ignored as they numerically and conceptually insignificant in 
this case. 

''The overall sign of h{z, s] is incorrect in [24]. 

®The incorrect argument in [27] stems from an incorrect relative sign between the "Inami-Lim" term and the QCD correction. 
^Cio is indeed very tiny in the SM as we shall show and is also stated in [27]. Hence, Cio(mc) = Cio(mT^). However we 
keep this contribution as it is important in Ap^ and has potentials of being largely enhanced by LHT as we saw in [18]. 
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2 SM Contributions to D ^ XJ+T 



SM contribution is dominated by Cg primarily with contributions from the purely electroweak part, coming 
almost entirely from the electromagnetic penguin, and an order of magnitude smaller contribution from the 
QCD correction coming from the four fermion operators Oi and 02- Contrary to what is stated in [27], 
Cj provides only a subdominant contribution in spite of its huge enhancement from the two loop 0{as) 
contributions. 

As we stated before, the SD contribution in SM is completely overshadowed by the LD contribution that 
comes from intermediate vector meson states. These resonance contributions lead to a branching fraction 
estimated in [24]: 

BRI^{D ^ Xue+e-) = BR^^\D X„e+e-) - 0(10^*^) (14) 

Rough estimates can show that the SD contributions for the branching ratio are much smaller than 
LD contributions. Why did we (and other authors) undertake a time consuming OPE analysis? Finding 
incorrect statements in published literature maybe be intellectually acceptable, even if such statements are 
of only academic significance. However we are driven by matters of much more practical interest: Some 
asymmetries on which distributions are based give a more direct access to SD dynamics, namely Apg and 
Aqp as we shall discuss next. 

2.2 Apg and A^p and Apg 

Asymmetries between — > X^l^l^ have a single source, namely AC — 1, AQ ~ 1 currents. Later we will 
comment on i?" — >■ Xul'^l^ vs. -> X^l^l^ , where — oscillations, in principle, could get involved. 
The normalized forward-backward asymmetry is defined from the double differential decay rate as 

, , /-I [-inno^ ^ Xjn-) £,TiD- ^ XJ+n] sgn{z)dz 

^FB S = i 15 

J_^[d^T{D± ^ XJ+l-)/dsdz]dz 
After performing the integral over the angular distribution we get 

-3 mcMC9i^^))s + 2^ic*Mc,{^l)] 



^fb(^) 



(1 + 2s) [iCifit + \Cwi^^)\ ) + 4 \Cj{nr (1 + f ) + {Cjifi)C^ifi)) 

(16) 



Since Cio is real, ^pb(^) picks up the real part of Cg and C7 which are both in general complex. Integrating 
over s we get 

A^B ^ 2 X 10-^ (17) 

j4pg(s) is mostly proportional to Cio which is tiny in the SM as expected from the suppression of FCNC 
in charm physics. Unlike the integrated decay rate, Apg is not very sensitive to infrared divergences. Since 
the angular distribution of the double differential decay rate is almost opaque to the LD contribution from 
SM^°, ylpg absorbs purely SM SD contributions - yet truly tiny! 
The CP asymmetry parameter Aqp{s) is defined as 

4rT(D+ XJ+l~) ~ 4rT(D- ^ Xul+l-) , , 

^ d_^riD+ ^ XJ+1-) + f.T{D- ^ XJ+1-) 

In general any Wilson coefficient Ci(/i) in the differential decay rate can be written as 

c^{^i) = + He 

where summation over j is implied. For T{D — > X^l^l^), Xj — > A*. The numerator will have contributions 
of the type 

^ 9(A'A'"*)3(e^4'"*) + 23(A')9(C°C-*) 
3? (C,(m)C*(/x)) ^ 3(A'A™*)3(C^er) + ^i^'M^lO + 3(A'"*)5(^°?f *) 



'A detailed argument on this can be found in [30] for the B mesons. A similar argument holds for the D mesons too. 
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2 SM Contributions to D ^ XJ+T 



Thc denominator will have contributions of the type 

3? ^ 3?(A'A"*)5R(^^er) + 3?(A')5R(ekf ) + 5R(A"*)5R(e°er) + (19) 

In the limit that only Cg has an imaginary component and terms proportional to V*ijV^f^ are ignored, the 
results stated in [31] are realized. The numerator is sensitive only to the imaginary contributions from the 
Wilson coefficients to the decay rate and the denominator to the real contributions as should be the case; 
for the relative phases in the matrix element is needed for a CP asymmetry. Unlike in the B mesons, in the 
case of 13 — > XJ'^l^ only Cio is purely real, and hence both Cr and Cg contribute to A^p. The integrated 
asymmetry in SM turns out to be 

_T{D+ ^ XJ+l~) -T{D- ^ Xul+l-) _4 

- T{D+ ^ XJ+1-) + r{D- ^ X.-J+1-) ^ 3 ^ 10 ; (20) 

i.e., still tiny. We also find that the bulk of A^^p comes from Cj and is due to the presence of the two loop 
0{as) contribution. Indeed, the Cg contribution, which comes from the mixing of Oi_6 with Og, serves only 
to suppress this contribution by an order of magnitude all of which stand in stark contrast to what happens 
in the analogous decay of the B mesons. Both of these contributions arise only after the inclusion of QCD 
radiative corrections and both of which are proportional to V*gV^g; ignoring the term proportional to V^*^!^^ 
as it is relatively much smaller. 

We can also look at Apg(s) which is the normalized difference in the forward-backward asymmetry in 
D XJ+l- and D -> Xal+l' defined as [32] 

Apq{s) = ^fb('^) + ^fb('^) ^21) 
^fb('^) ~ ^fb('^) 

In the limit of CP symmetry ^pg(.s) and Ap-g{s) have to be exactly equal in magnitude but with an opposite 
sign [31,33]. As the forward-backward asymmetry is defined in terms of the positive anti-lepton, Ap^{s) and 
j4pg(s) have opposite signs, ylpg(s) is sensitive to the phase in C7, Cg and Cio- The SM offers phases only 
in C7 and Cg in Z) — >■ X„Z"*"Z~ and none in Ciq. Hence the integrated asymmetry turns out to be tiny. 

A^^{s)ds = ^ 3 X 10-^ (22) 

If NP brings about any new phases in either C7, Cg and Cio, ^fb stands a chance of large enhancements. 

Putting together Apg , A^p and Apg gives us a good insight into the sizes of the phases in the Wilson 
coefficients, ^p^ is sensitive to the size of Cio and the real parts of C7 and Cg, while A'^p gives us an idea 
of the size of the phases in C7 and Cg and Ap^ is sensitive to phases in all C7, Cg and Cio- Within the SM, 
we can conclude from our numbers, the size of Cio is extremely small, as it should be, since it suffers from 
a very strong GIM suppression, and it also lacks a phase. Both C7 and Cg have phases because of the QCD 
corrections which are more prominent in the case of D mesons than in B mesons as the purely electroweak 
contribution is truly tiny. Along with the decay rate, these provide us a very useful tool to probe into the 
flavour structure of any NP models and new sources of FCNC. 



2.3 A Note on LD Dynamics' Impact on the Asymmetries 

The asymmetries discussed in the previous section do not incorporate any SM LD contributions, neither 
in their extraction from the difference in the partial decay rates nor in their normalization. This might be 
unreasonable considering that SM SD contributions yield a branching ratio of only 1.5 x 10~^, while SM LD 
yields something like O(10~^) even considering that LD estimates come with very large uncertainties. Let 
us add a comment why LD physics have only a little impact on these asymmetries, either from the difference 
in the numerator or the normalization in the denominator. 

For A'^p and Apg it is obvious that SM LD physics cannot contribute a CP violating phase, and any CP 
asymmetry in this process has to come from SD physics whatever the origin is - SM or NP. Ap^ is more 
sensitive to SM LD 'pollution' in its definition as the difference in the hemispherical integral of the double 
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3 On LHT Scenarios 



differential decay rate. SM LD contributions to it will primarily come through final state interactions and the 
dominance of light internal quarks in this process and quantitative statements on their sizes require separate 
analyses on each exclusive process. SM LD physics can also make their presence felt in all these observables 
through the normalizations entering the definitions of these asymmetries. A natural way to remove such LD 
contributions is to cut off the dilcptonic mass distribution around the p, w and widths. 

Making such cuts we find that it will decrease A^^ and increase A^p and by around 10% to 20%. 
Not surprisingly, these cuts affect Ap^ more than A^p and Apg. Similar effects arise, when we include LHT 
contributions as discussed below. More sophisticated cuts could yield even better results. As pointed out 
below, only Apg has a realistic chance to be measurable by LHCb and a Super-Flavour Factory. It seems to 
us that such a theoretical uncertainty is fully acceptable for a search for NP as of now. 

As a final point: including non-rcsonant LD effects and making a cut on the resonances seem to produce 
opposing effects and can well nullify each other; however only a detailed study can resolve this issue. 

Throughout the rest of the article we will stick to the definition of the asymmetries in terms of SD 
operators and infer on our results accordingly keeping in mind that a sufficiently motivated reader will 
already have had gone through this section by then. Any NP contribution to these asymmetries are strictly 
SD. 

2.4 Comments on D°/D° XJ+l' 

The branching ratios for neutral D — >■ Xul'^l^ are again dominated by SM LD contributions and their size 
is comparable to what is stated in Eq.l4, namely of order 10~^. For Ap^ and A'^p one has a much more 
complex system in hand for neutral D mesons, because D'^ — oscillations have been found on the level of 
0.5% — 1 % for XD and j/d, see Eq.l, which might be still consistent with the SM. Since the SM asymmetries 
Ap-Q and A^p are so tiny, the — oscillation 'background' is irrelevant. 

3 On LHT Scenarios 
3.1 The Flavour in LHT 

The SM predictions presented above leave a large range in rates for these rare transitions, where NP could 
a priori make its presence felt. So-called Little Higgs models mentioned in the Introduction have been 
studied extensively over the past decade as a possible NP scenario [8,9]. There the Higgs boson appears 
as a pseudo-Nambu-Goldstone boson of a spontaneously broken global symmetry. Rather than attempting 
to solve the hierarchy problem, they 'delay the day of reckoning' and address a maybe secondary, yet very 
relevant problem, namely to reconcile the fact that the measured values of the electroweak parameters show 
no impact from NP even on the level of quantum corrections with the expectation that NP quanta exist 
with masses around the ITeV scale so that they could be produced at the LHC. In some of these models, to 
achieve this program consistently, one needs an additional discrete symmetry [13-15], called T parity. One 
way of consistently implementing T parity also requires the introduction of the so-called 'mirror' fermions - 
one for each quark (and lepton) species - that are odd under T parity and familywise mass degenerate. This 
introduces two 3x3 mixing matrices Vud and Vhu neither of which need to be close of the CKM matrix, 
but they are related to each other [34]: 

yLyHu = vcKM (23) 

Since the CKM matrix does not differ too much from the identity matrix, one realizes that LHT contributions 
exhibit a clear correlation of the phases in the charm and strange sector. 

In this note we will analyze a subclass of Little Higgs models, namely Littlest Higgs Models with T parity 
(LHT) [11. 16]. In our view they possess several significant strong points: 

• They contain several states with masses that can be below 1 TeV; i.e. those states should be produced 
and observed at the LHC. 

• Compared to SUSY models they introduce many fewer new entities and observable parameters. 

• Their motivation as sketched above lies outside of flavour dynamics. Thus they have not been 'cooked 
up' to induce striking effects in the decays of hadrons with strangeness, charm or beauty. 
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4 Numerical Findings on LHT Contributions 



• Nevertheless they are not of the minimal flavour violating variety! 

• The impact of LHT dynamics on K, B and also D transitions has been explored in considerable detail, 
and potentially sizable effects have been identified [18,35-38]. 

• Especially relevant for our study is the fact that they can have an observable impact on D° — _D° 
oscillations [17,20]. Also sizable indirect CP violation can arise in decays [17] very close to the 
present experimental upper bounds. Having seen such large effects in AC = 2 transitions coupled 
with the possibility of the existence of large CP violating phases, one would naturally ask whether it 
is possible to see the same in AC = 1 transitions such as Z? — >■ 



3.2 LHT contributions to D ^ XJ+l' 

As T parity forbids tree level coupling of SM particles with the new T-odd particles, LHT makes it presence 
felt only through loop contributions from internal mirror fermions and heavy gauge bosons. Unlike in the 
case of B and K, the new T-even heavy top quark does not contribute and hence any new contribution from 
LHT comes from the T-odd particles. The following are the modifications of the SM functions in Eq.5. The 
auxiliary functions are defined explicitly in Appendix B. 

Ci{x) = (^xS{x) ~ 8xR2ix) + ^x + 2xF2{x) 



^i(^) = \ j2 {d^{x)+Qu\eo{x)+Qu^E^{x)^ 



,2 



^i(^) = \j2 (d',{x)+Q^\e',{x)+Q^^E'^{x')^ 
E[{x) = ~ [e',{x) + \e',{x) + l^E',{x')^ 



Y,{x,y) = ^^[5(a;)+Fv^(x,y)-4(Gz(x,y) + G^(x',y') + G,(x,y))] 



(24) 



where 
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tan^ 
1 



y' = ay , 7] ^ - 
a 

(25) 

= 2/3 is the charge of the up-type quarks, m//^ is the mass of the mirror quark in the i*'* family and 
m;^ is the mass of the heavy internal neutrino. The functions F\^{x,y), Gz{x,y), GA{x,y) and Gri{x,y) 
are contributions from WW, ZZ, AA and ZA box diagrams with heavy internal neutrinos respectively. 
A complete list of Feynman diagrams can be found in [38]. Since the operator structure is the same in 
LHT as in SM, the expression for the decay rate and the asymmetries remain the same with the necessary 
modifications of the Wilson coefficients. QCD corrections to the LHT contributions have been ignored; after 
all we do not know the model parameters, and these numerical exercises serve to show whether such models 
can be significant for such observables. 



4 Numerical Findings on LHT Contributions 

The structure of the mirror fermion sector leaves us with a lot of liberty to choose the parameter space we 
wish to scan. However, constraints from B and K physics set very stringent limits on the viable parameter 
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space for probing D physics. In what foUows, we define an operational parameter space and what effects 
LHT can bring about in _D ^ Xul^l^ ■ 

4.1 LHT Parameter Space 




mHi(GeV) 1000 mHi(GeV) 



Figure 1: Parameter space of the mass of the mirror quarks 




in our previous work [18]. We vary the mirror fermion masses and the mixing angles and phases over 
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parameter sets keeping the breaking scale of the non-linear sigma model fixed at 1 TeV^^. The LHT has 20 
new parameters of which the ones which will be relevant to us are as follows: 

• The LHT breaking scale / = 1 TeV is fixed by choice. 

• The masses of the three familywise degenerate T-odd mirror quarks, m/fi, 'T^^ffa range from 300 
to 1000 GeV. 

• There arc three independent mixing angles in Vhu, ^121^137^23- 

• There arc three irreducible phases in Vhu-, ^121 ^131 ^23' 

The parameter space used for these analyses is a set that satisfies all experimental constrains from B and 
K physics. A small parameter set was also used which did not follow such constraints to check whether 
constraints from B and K physics affects LHT contributions to D physics. However, even the parameter set 
that is not constrained does not have large mass hierarchies in the mirror fermion sector. 

The mass spectrum for both the parameter sets is illustrated in Figs.l. Using Eq(23), the angles and 
phases of Vhu were calculated from those of Vnd and hence were constrained by B and K physics too for the 
first parameter set and not so for the second. Histograms of the parameter space of the angles and phases 
are shown in Figs. 2. The angles and phases are family-wise paired. 

4.2 Impact on TgD (D XJ+l') 



3500 p 
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(a) Constrained by B and K Physics 



m,„= 400 GeV 




m,„= 600 GeV 




m(„= 1100 GeV 



(r,ot-rsM)/rsM % (SD) 



(r,<„-rsM)/rsM % (SD) 

(b) Not Constrained by B and K Physics 



(r,ot-rsM)/rsM % (SD) 



Figure 3: Percentage change of the SM SD contribution to TsriiD Xul^l ) due to LHT effects. 

As we have seen in Sec. 2.1, the dominant SM contributions to Fsd {D — > Xul^l~) are through the 7 
penguin in the Og operator and some from the mixing of (3i_6 with Og. The subdominant contributions 
come from the two loop 0{as) term in Cj. Hence part of the dominant effect and the subdominant effect 
both come from QCD corrections. All other contributions are smaller by orders of magnitude. We saw in 
our previous work on D 77 and D /^^M^ [18] that LHT is capable of producing large enhancements 

^ ^ An analysis of the parameter space can be found in [39] . 



12 



4 Numerical Findings on LHT Contributions 



through box diagrams with internal heavy fermions and heavy gauge bosons just as in D^-D° oscillation [17] 
and somewhat moderate enhancements to penguins. However, the enhancement to effective 7 vertices 
are tiny compared to the SM contributions. 

In D — > Xul^l^ we confirm our previous conclusions. The decay rate was calculated for three different 
internal heavy neutrino mass of 400 GeV, 600 GcV and 1100 GeV. We see very tiny change in both the 
differential decay rate and the integrated decay rate at 0(1%) as can be seen from Fig. 3 in which the abscissa 
represents the percentage enhancement to the SD decay rate after inclusion of LHT. Removing constraints 
from B and K physics does not make much difference either as can be understood from comparing Fig. 3 (a) 
and Fig. 3(b). Our result is different from what was found in [40] which used the Littlest Higgs model 
withoMt T parity. Without T parity the SU{2) custodial symmetry protecting the electroweak p parameter 
is explicitly broken at scales below 4 TeV [12] which is brought about by the U{1)h gauge boson. This 
model allows large tree level FCNC mediated by the coupling of the Zh and Ah heavy gauge bosons with 
the SM quarks. The enhancements seen in [40] is because of the existence of these tree level FCNCs which 
are absent from LHT. Hence, LHT makes almost no changes to C7 or Cg above SM contributions and hence 
fails to enhance the decay rate. This essentially means that any LHT contribution to the decay rate fails to 
significantly dent SM SD contributions and is completely swamped by SM LD effects. 



4.3 Impact on A 
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Figure 4: Ap^ after including LHT effects. 



As we have seen earlier in Sec. 2. 2. the SM contribution to A^^ is all but nonexistent due to tiny SM 
contributions to Cio . In LHT Cio gets enhanced by orders of magnitude which brings about a large enhance- 
ment in ApQ. This effect is similar to what was observed in [40] but comes from box diagrams involving 
T-odd heavy internal degrees of freedom rather than tree level FCNC. It is also commensurate with the 
enhancement we found in SD contribution to T{D° — >■ fi^ [18] from LHT. This is due to the fact that 
SD contribution to r(D° — > jX^ p~) comes from Oio and Apg is highly sensitive to the same. However, even 
with such a large enhancement, the absolute value A^^ after including the LHT enhancement can at most 
be of 0(0.5%) as can be seen from Fig. 4. Studying Fig. 4(a) and Fig. 4(b) shows that removing constraints 
from B and K physics creates large relative enhancements - it can be as much as 1% (or rarely more) - but 
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(b) Not Constrained by B and K Physics 

Figure 5: Enhancement to over SM after including LHT effects. 
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Figure 6: Large enhancements in after including LHT effects 
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they do not enhance it to sizable absolute effects. 

On the other hand A^p depends mostly on C^ and Cg, which we have already seen suffers almost 
no enhancement from LHT. However, as A^p is quite sensitive to the phases in these coefficients, it gets 
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enhanced by a few factors over the SM value and can be as large as four times the SM value. This can be 
seen from Fig. 5 where wc plot the ratio of the total A^p including the LHT enhancements to the SM value 
of the same. However, this still keeps A^p at O(10~*) — 0(10""^) and hence the absolute measure of the CP 
asymmetry is still experimentally challenging. The unconstrained parameter set allows for slightly larger 
enhancements to but is limited to almost the same order of magnitude. 

The contributions from LHT models can enhance Apg so much as to bring it up to possibly measurable 
values. As pointed out in Sect. 2. 2, Ap^ is sensitive to any phase in Ciq. LHT can not only enhance the 
magnitude of Cio by orders of magnitude, but also brings about the possibility of existence of a very large 
phase in it. For Cr and Cg the effect is dominated by phases from QCD radiative corrections. The existence 
of this large phase and the tangential dependence of Ap^ on it results in the huge enhancement that we 
see in Ap^ as illustrated in Fig. 6 where we plot the total Apg after the inclusion of LHT effects. This 
is commensurate with what was observed in [17] for CP violation in D° — I)° oscillations. For both the 
constrained (Fig. 6(a)) and unconstrained (Fig. 6(b)) sets more than 10% of the parameter set can produce 
asymmetries of 0(10%) or greater! As explained in Sect. 2. 2 we relate Ap^ with SD contribution to — 
l~^l^X, which amounts to a branching ratio of 1.5 x 10~^. With a sample of 10^"^ D mesons NP intervention 
should be measurable for Apg > 5%. Such effects could be within the reach of to the currently running 
LHCb experiment and ones like the planned and approved SuperB Collaboration. 

4.4 Correlation between Apg and Fsd {D — >■ Xul~^l~) 
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Figure 7: Correlation between Ap^ and Fsd {D — )• X^l^l ). The gray points represent the constrained set 
and the black ones for the unconstrained one. The big black spot represents the SM values. 

In Fig. 7 we plot the correlation between Apg and Fsd [D Xul^l~). These plots look at first chaotic, 
yet a careful (and time consuming) study reveals a pattern. For 'low' heavy neutrino masses mi^ ~ 400 
GeV, 600 GeV the LHT parameter sets that enhance Fsd {D — > X^l^l^) can produce positive Ap^ in some 
regions and negative ones in others; on the other hand sets decreasing Fsd {D — > Xul^l') can also produce 
positive and negative ^p^, but in others regions. Those 'low' heavy neutrino masses are within the range 
of masses used for the mirror quarks in this study. However, for nii^ = 1100 GeV those LHT parameters 
that increase the SD branching ratio produce mostly a negative Ap^ , while sets decreasing the SD branching 
ratio lead mostly a positive Ap^. At this mass the mirror neutrino is heavier than any of the mirror quarks. 

This behavior can be understood quite well. Dependence on the heavy neutrino mass exists only in the 
box diagrams. The sign of Ap^ depends on the sign of Cio which depends on the relative size of the mirror 
quarks and heavy neutrino masses which comes from box diagrams. At less than 1 TeV, this can go either 
ways with the neutrino being either heavier or lighter than one or more of the mirror quarks with slightly 
greater chances of being lighter than them. At above 1 TeV the heavy neutrino is always heavier than 
the mirror quarks used. Also enhancements to Fsd {D — >■ Xul~^l~) are quite sensitive to the box diagrams 
through both Cio and Cg. This leads to the sharp change in the correlation we see in Fig. 7. 

Wc do not see any such correlation in vs. Fsd {D Xul'^l') or Ap^ vs. A^p; after all A^p is blind 
to Cio and hence less sensitive to the box diagrams as in this case they manifest themselves only through 
Cg. Furthermore we see no correlation in Ap^ vs. Fsd {D — >■ X^l^l^) or Apg vs. Ap-Q since Ap^ is not 
affected by the magnitude of Cio, but rather by the CP violating phase in it; that information is lost in both 
TsD {D — >■ Xul^l~) and Apg. Lastly, we see no correlation between ^^p and App even though both are CP 
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violating parameters. This reinforces our earlier statement that the sources of CP violation are distinct in 
these parameters with the former coming from phases in Cy and Cg and the latter coming from a phase in 
Cio- 

5 Further Insights into FCNCs in LHT-hke Models 

Our work on the impact of LHT on D'^ — > fi'^ fi~ [18] had lead us to some general conclusions on the 
structure of FCNCs within a LHT-like framework. As defined previously, this framework contains 

• A second sector of fermions that arc an exact copy of the SM ones. 

• Mass mixing matrices which are unitary and loosely connected to Vckm (Eq.23). 

• Possible large angles and phases in the mass mixing matrices. 

• Possible large hierarchies in the masses of the mirror quarks. 

• A symmetry, like T parity, segregating the NP sector from the SM sector, hence forbidding tree level 
FCNC. 

We shall, after further investigation, relax the second condition to: 

• Mass mixing matrices that are constrained by a relationship between the one(s) connecting the new 
Up-type quarks with the SM down-type quarks to the one(s) connecting the new Down-type quarks 
with the SM up- type quarks. 

In general, FCNCs are a very sensitive probe to the details of the flavour structure of both the SM and 
any NP models as they highlight not only mass hierarchies within a theory but also are sensitive to phases 
within the same. Moreover, it is possible to disentangle the effect of phases and fermionic mass hierarchies 
on FCNCs in a model independent way if we have access to more observables. 

D — >■ Xul^l~ wins over — > 77 and l^^f^' by leaps and bounds in this respect. While — >■ 77 

is sensitive mostly only to O7 and — > fi^fi^ is sensitive only to Oio, D — > Xul~^l~ is not only sensitive 
to all of that but also to many more. Moreover, the final state being a three body final state, this channel 
can also be probed through forward-backward and CP asymmetries hence opening the possibility of probing 
phases in any model, too. 

As we noted above, the SD contribution to branching fractions is dominated by the photonic penguin in 
Og while Apg is highly sensitive to Oio, A'^p to the mixing between O7 and Og and Ap^ to the phase in 
CiQ. Studying D — > Xul~^l~ in a sufhciently precise way, we can learn the impact of several operators and 
then comment on other rare decays. Logically we should have started our analysis with D — > Xyl^l^ and 
then applied our findings to the simpler cases of two-body rare decays. Instead we started with our analysis 
of two-body rare decays, from which we extracted some conjectures; they happened to be correct in more 
general theoretical situations. 

From our results wc sec that, through LHT dynamics, Apg gets orders of magnitude enhancement through 
the enhancement of Cio which is commensurate with the orders of magnitude enhancement that we had seen 
in the SD contribution to which was ultimately overshadowed by the LD contribution to the 

branching fraction. We see almost no enhancement to C7 which agrees with the lack of enhancement that 
we noted in D'~^ — > 77. It is not that C7 gets absolutely no enhancement from LHT. The purely electroweak 
part of C7 does get moderately enhanced. However, this enhancement is completely overshadowed by the 
SM two-loop 0(q!s) QCD correction, something which is peculiar to the D meson system and not seen in 
the B system. In addition, we also see almost no enhancement to Cg coming from the fact that photonic 
penguins are not enhanced by LHT, which results in a lack of enhancement to the branching fraction of 
D XJ+1-. 

Intuitively one might think that the situation for A^p should be different as it is sensitive to phases in 
Cg and C7, and LHT allows for large phases in the mixing matrices. In a purely electroweak SM scenario 
this would have been true. However, gets an unusual boost within the SM from the unusual two-loop 
0{as) QCD correction, something that LHT can barely overcome. Coupled to the fact that Cg does not 
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gain much from LHT, enhancements to A^p fail to impress. The validity of the previous statement is further 
tested when we see orders of magnitude enhancement in Ap^ as it gets a boost from the introduction of a 
large phase from LHT while it is completely clean of phases from the SM. 

So what does this tell us about LHT-like models and their effects in flavour physics? In principle in 
any such model, large angles and phases in mass mixing matrices and large fermionic masses and sharp 
hierarchies amongst them are possible. However, one has to satisfy the experimental constraints that we 
already have in B and K physics. At this point one will have to choose between large angles and phases or 
huge hierarchies in the fermionic masses as experimental data already tell us that the extra fermions have 
to be heavy. Making such a choice automatically limits the size of NP intervention in yet unobserved FCNC 
processes, specially in AF = 1 processes even if AF = 2 processes can escape these limits and absorb NP 
contributions. However, an exception to this rule occurs when these large phases from New Physics are laid 
bare and have purely SM electroweak effects to compete with as in the case of Apg where we sec large effects 
even in a parameter which is a measure of a AF = 1 process. Hence, let us have a more thorough look at 
the diagrammatic details of NP intervention from LHT-like models. 

6 On Boxes and Penguins in New Physics^^ 

Due to CPT symmetry, CP violation can enter only through complex effective couplings. In the SM they 
can arise only for the weak boson couplings to quarks as described by CKM matrices. They are necessarily 
unitary, since all quark masses are given by a single VEV of a neutral Higgs multiplied by numbers, not a 
matrix. The concept of 'weak universality' was first put forward by Cabibbo in 1967 [41]. Afterwards it 
was scrutinized experimentally. Later it was understood if the weak forces are embedded - as it applied to 
the SM - in a single non-abelian gauge theory, weak universality has to hold. There is the Singular Value 
Decomposition theorem which tell us that the matrices relating mass and flavour left Up and Down quarks 
are unitary and therefore their matrix product - the CKM matrix ~ is also unitary. For N = 2 families 
their phases can be transformed away, for iV = 3 there is one irreducible phase that is therefore observable. 
CP violation can surface in processes where quarks from three families can contribute as real or virtual 
entities. The latter happens in SM, since FCNCs arise effectively through quantum corrections, namely box 
and penguin diagrams. There are box diagrams in only one kind, the WW box. Penguins come in three 
varieties: the Z^^ the 7 and the gluon penguins are particularly essential for direct CP violation. 
For N families the unitarity of Vckm translates into^'^ 



with |Ai| > for i = 1, 2...N families with {i} coming from the down-type sector and {j, k} coming from the 
up- type sector - i.e., a triangle equation in the complex plane: CP violation arises at the mass generation. 
If two of the quark were mass degenerate, CP invariance would survive mass generation - yet it is not the 
case on our world. 

If the SM contained a fourth family, it could - and probably would - have two more observable phases 
with more independent CP asymmetries. 

The very fact that in the SM one has to study transitions where quarks from three families contribute 
reduces the number of 'interesting' cases for CP violation very significantly. The observation that absolute 
values of off-diagonal elements of Vckm are small, even tiny, in most such cases are experimentally 'chal- 
lenging' at least - except in beauty decays, since leading decays are so suppressed and therefore make them 
experimentally challenging for a different reason. It was expected that absolute values of the off-diagonal 
elements are small, but it was surprising that |Vf,c|, \Vbu\ etc. are so tiny, namely much smaller than \Vsu\- 
Originally it was conjectured - based on no good theoretical reason - that also the CKM phase is small. 
Now we know that the CKM phase is not small. Therefore NP is likely to exhibit also sizable phases - like 



^^and Seagulls too, but for now we shall ignore them. 

^^This triangle relation holds when {j, k} belong to the up-type sector. If they belong to the down-type sector we have 
Xi = V^j* V^j, instead with {«} coming from the up-type sector. All other arguments hold true. 



N 




(26) 



i=l 



LHT. 
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The size of CP asymmetries depends on the phases, the absolute values of the quark mixing matrix 
elements and the masses of the internal virtual gauge bosons and fermions and their mass hierarchies. Let us 
analyze to what degree each criterion can be fulfilled in this class of NP models we call LHT-like for box and 
Penguins. The range of masses of mirror quarks and Wh and Zh is rather limited - for two very different 
reasons: (a) None have been found up to ~ 200 GeV. (b) They could be much higher like 10 TeV or even 
100 TeV. Then they could not be directly produced at the LHC. Therefore we stop our analysis at the 1-2 
TeV scale. In that case some of the new states could mix significantly. The gauge boson mass is also of 
the 0(1 TeV). Large angles and phases are possible but limited by experimental constraints from B and K 
physics. Also, quite important to our analysis, hierarchies in fermionic masses are very small, a lot smaller 
that what we find in the SM. 

In LHT-like models the operators with new degrees of freedom are similar to the SM SD ones. Hence 
they scale similarly to the SM SD operators which has been established for decades now [22] and has been 
used extensively to qualitatively judge the size of Savour dynamics ever since. Let us reexamine this scaling 
behavior. Of course, this scaling is only an approximation of the detailed formfactors, however, most of 
the time such approximations arc enough to estimate the size of many effects in fiavour physics. The mass 
hierarchy in the new fermionic sector can be established with 

mf = mj/if yi = l...N (27) 

Here a distinguishes amongst the members of each family which defines the the sector to which the ferniion 
belongs and it is not summed over. It is possible to set ^ 1 and choose to" to be any finite mass 
representative of the phyics in consideration. The formfactors in flavour physics can be approximately 
expressed as 



F{x) = /„(x)(log(x))" where, /„(x) = x", ti e Z, m = 0, 1 (28) 

Ignoring QCD (or QCD-like) corrections, any matrix elements involving the processes mentioned above are 
of the form 

N 

A^~^A,F(x,) (29) 

i=l 

Here, xf = {ml /mcY is the commonly used square of the ratio of the internal fermion mass to the mass scale 
of the massive gauge bosons. The superscript a has been dropped as the matrix element usually involves 
just one sector. Using the unitarity relation in Eq.26, it can be shown that under the hierarchy defined in 
Eq.27, Eq.29 scales as 

N N N 

KF{x,) = /„(xi) ^ KF{hl) + F{x,)Y^ Kfnih^) (30) 

i—l i—1 i—1 

Note that the mass scale of the fermions docs not suffer from unitarity suppression but the hierarchy does. 
In other words, the only two ways of getting large matrix elements are either to start at a very high mass 
scale for fermions or to build a very strong hierarchy that will illudc the unitarity suppression, or both. 
Now, let us have a look at boxes and penguins. 

Box diagrams 

Box diagrams scale as x for both large and small x. Setting n = 1 and m = in Eq.28 we see that 
contributions from these will scale as 

JV 

X-a;i^A,/i2 (31) 

4=1 

This clearly tells us that box diagrams are sensitive to both mass scales and large hierarchies. Hence 
any NP model containing either fermions with masses comparable or greater than the gauge bosons 
or having a large hierarchy amongst the families or both will make large contributions through box 
diagrams. 
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Zh Penguins 

Z]^ Penguins scale as x log(2;) for small x and as x for large x. In the regime in which they scale as x 
the conclusion is the same as above. Setting n = 1 and m = 1 in Eq.28; we see that if, for some reason, 
NP offers small x, the contribution will scale as 



Hence, in addition to the conclusions drawn for large x, Zj^ Penguins will get large contributions from 
NP if the fermionic masses are light compared to the gauge boson mass and shows large hierarchies. 
However, such scenarios will be very unusual if not unheard of in the NP models currently under 
consideration. 

7 and chromomagnetic Penguins 

Electromagnetic and chromomagnetic penguins scale as x for small values of x but are asymptotic to 
a constant for large values of x. Hence they can only benefit from large hierarchies at fermionic mass 
scale smaller than the gauge boson. Even if NP has large fermionic masses, even with large hierarchies 
to offer, it will show up only as moderate enhancements in this class of Penguins. 

On the other hand, photonic Penguins scale as log(.T). Hence, n = and m = 1 in Eq.28, we see its 
contribution to matrix elements will scale as 

N 



NP intervention will fail to produce any enhancement in photonic Penguins unless it has very sharp 
hierarchies to offer in its fermionic masses. Large mass scales will have no effect on photonic Penguins. 

A look at the mass spectrum of the heavy fermions, Fig.l, which we used in the previous study [18] 
and in the present one, show that we have heavy mass scales in the spectrum, but not large hierarchies. 
As a result, we get large enhancements to purely electroweak processes which involve box diagrams and 
Zh penguins. Such is the case for — > and Apg where we see orders of magnitude enhancement 

to SM SD rates and also in — oscillations. Apg too benefits from this effect. For electromagnetic 
penguins and chromomagnetic penguins, NP intervention can only be moderate. This explains the small 
enhancements that we saw in D° — > 77 which is primarily sensitive to O7. However, in Z) — > Xul~^l~, the 
dominant contribution, by orders of magnitude is the photonic penguins in O9 if we ignore QCD corrections 
which are anyways blind to NP. Even if LHT manages to enhance the other contributions, it falls short of 
the SM contribution to Cg . Even the LHT contributions to the box diagrams and Zl Penguins in Og fails to 
overcome the SM contribution to the photonic Penguins. Moreover, the photonic Penguins do not see much 
NP intervention as we have already argued. 

This gives us a qualitative way of understanding why LHT fails to enhance processes driven by AC = 1 
dynamics involving the D meson system unless QCD effects are completely absent. This argument can 
potentially be extended to any other model which shares the same flavour structure as LHT, what we have 
previously defined as LHT-like. Of course, we always have to keep in mind that we are heavily constrained 
from B and K physics and that limits us in the parameter space that we can choose to work with. 

7 Conclusions 

Charm hadrons stable under strong forces were predicted to keep the SM consistent with the observed 
suppression of strangeness changing neutral currents, as were their preference to decay into strange hadrons 
with decays into non-strange hadrons being Cabibbo suppressed and to preserve renormalizability; they were 
even found in the expected mass range. Therefore decays of charm hadrons were hardly seen as worth probing 
for manifestations of NP. Furthermore one realized that kaon and even more, beauty hadrons could clearly 
exhibit NP signals, since their leading SM decays are Cabibbo and KM suppressed. A minority of authors 
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argue that the probe for NP should not be given up in charm decays, since the SM weak phenomenology is 
'dull' and one need 'only' much more statistics there. 

The case was somewhat strengthened after the observation of D° — oscillations, although it is not 
outside some SM estimates. Furthermore it was found that a class of NP models like LHT that is motivated 
from outside the flavour dynamics, can yet produce a sizable contribution to the observed D'^ — D'^ signal 
and can create much stronger indirect CP violation [17]. 

D Xul^l^ is a good laboratory for NP as it is rich in its operator structure and involves almost 
anything that FCNCs have to offer, including direct CP violation. However, we find that LHT fails to create 
a significant dent to SM SD contributions to this channel other than significantly enhancing Apg and Apg 
for reasons we have explained above. Moreover, SM LD contribution dominates over both SM SD and LHT 
contributions to the branching fractions. We reiterate the conclusion that we had come to in our previous 
work [18]: while LHT can contribute significantly to AC = 2 processes, it fails to dominate in AC = 1 
processes with orders of magnitude enhancements unless the contribution appears through box diagrams 
and are bereft of relatively large QCD effects. 

We also go ahead and have a second look at what we previously defined as LHT-like models: again 
we found that certain conclusions can be drawn about NP's flavour structure by studying weak decays 
independent of the other details of the underlying model. Absence of large hierarchies and unitarity of the 
new mass mixing matrices within such NP models heavily limits the new FCNCs. Also, we have shown that 
experimental limits in B and K physics can be directly ported to give constraints on NP intervention in D 
physics. 

During our analysis we developed a more general conjecture: if NP models affect the dynamics of both 
the up- type and down-type quarks in a tightly correlated way - as shown in Eq.(23) in the case of LHT - 
they will mostly fail to contribute significantly in AC = 1 dynamics even if they play a major role in AC = 2 
processes. For constraints from AS" = 1 and A_B ~ 1 reactions will suppress AC = 1 coupling greatly, since 
sensitivity for NP is often greater in K and B decays because of their leading SM transitions arc Cabibbo 
or KM suppressed; otherwise it would oversaturatc for AC = 2 effects. The crucial feature on the right 
hand side of Eq.(23) is only that the CKM matrix is very close to the identity matrix; the same conclusion 
should apply for any matrix in the right hand side that is close to the identity matrix. Of course, the validity 
of this conjecture is yet to be tested in other AC = 1 processes like direct CP violations in nonlcptonic 
charm decays, keeping in mind that charm transitions could still produce surprises for us - and about SM's 
'ability' to cope with them. In the event that future experiments reveal a clear manifestation of NP through 
enhancements in the AC = 1 processes, it is unlikely that any LHT-like model can be a 'culprit'. 



8 "Postmortem" 

People reading this paper closely will realize how much theoretical working was needed, yet hardly any useful 
results were found at the experimentally observable level. They will ask, "what did you learn from your 
pain?" . One author, who did most of the work, probably asks himself the same question. One of the other 
authors will reply with a typical German answer: "You learn so much about NP models and field theory." 
The "long suffering" author thought of Robert Herrick, a disciple of Ben Jonson: 

NO PAINS, NO GAINS 

If little labour, little are our gains: 
Man's fortunes are according to his pains. 

- Hespendes 752, (1648) 
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APPENDIX 



Appendix A: QCD Corrections 



Here we include some of the numbers and functions which appear in the QCD correction of the operators. 
The function f{z) in Eq.A.l is that which appears in the 2-loop 0{as) correction to C7 is given by [23]. 



1 

"243 

+108 \og{zf + 36 log(z)3 ) z + ( 324 - STOtt^ + (1728 - 2167r2) log( 



576n^zi 



+324 log(z)^ + 36 log(z)-' F + 1296 - 127r^ + 1776 log(z) - 2052 log(z) 



3672 - 2887r2 - 1296C(3) + (1944 - 3247r2) log(z) 

2 I /'I'TOO oi(;_2\ 



(lU - 6tt^ + 18 log(z) + 18 log(z)2^ 2 + - 54 - 67r2 + 108 log(2) 



+18 log(2)2) z^ + f 116 - 96 log(2)) 2^1 + o{z^) 



(A.l) 

The vectors a and z [23] which appear in C7 are given by Eq.A.2 and Eq.A.3. The matrix X is given by 
Eq.A.4. 



X = 



M' S' ^' ~i' -0-4230, -0.8994, 0.1456 

M' S' ^' -0.8451, 0.1317 

/ -3.5687 2.5813 0.4 0. 0.6524 -0.0532 -0.0034 -0.0084 \ 

-4.0742 2.7827 0.4 0. 0.8461 0.0444 0.0068 -0.0059 

-22.423 18.290 0. 0. 4.3019 -0.1241 0.0001 -0.0452 

-23.434 18.693 0. 0. 4.6894 0.071 0.0206 -0.0402 

9.8081 -8.8366 0. 0. -0.7779 0.0289 -0.0486 -0.1739 

\ 2.8271 -3.2361 0. 0. 0.4903 0.0433 -0.1303 0.0056 / 

The Wilson coefficients of the operators Oi^g at /i = nib are given by [42] 

8 



with the initial condition 



(A.2) 
(A.3) 

(A.4) 



The coefficient matrix hji given by 





i=3 


-2)<- Vj = 1 . 


..6 




C{mw) = 


{0,1,0,0,0,0} 




given by 








/ 1/2 - 


-1/2 







1/2 


1/2 








1/6 - 


-1/14 0.0510 


-0.1403 -0.0113 


0.0054 


-1/6 - 


-1/14 0.0984 


0.1214 0.0156 


0.0026 





-0.0397 


0.0117 -0.0025 


0.0304 


V 


0.0335 


0.0239 -0.0462 


-0.0112 / 



(A.5) 



(A.6) 



(A.7) 



The Wilson coefficients given in Eq.A.5 can be calculated using the anomalous dimension matrix given 
in Ref. [23] with five active flavors for /i > m^. To calculate C(7tIc) the same procedure with four active 
flavors for > /i > rric will ycild the necessary results. 



Appendix B: Auxiliary Functions 

The auxiliary functions that come into the LHT contributions to — > Xul~^l~ are listed below. They can 
also be found in [38]. 
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Gz{x, y) = -|;7(x, y, 1) Ga{x, y) = ^Gz{x, y) 

Gr,{x,y) = --^U{x,y,i]a) Fw{x,y) = ^x - Fr,{x,y) + 7F(i{x,y) + m{x,y,l) 

J'2[X) ~ -2\ {1-xyi + ~J ^5[X,y) - - 

C/'^-i — ^ ( x^-2x+4 \ , 7-3: A TT( \ _ x^ log(x) log(i/) t ^osiv) 

b(X)-X\^ (1-^):. iOgi^J + 2(Ti:jy j U[X,y,Tj) - + (l__.y)(^_.y)(;,_:r) + (l-„)(3:-„)(„-y) 

(B.l) 

The function iS'(x) is a contribution from the penguin diagrams with internal mirror quarks which 
was pointed out in [43] and subsequently in [44]. It replaces the divergence mentioned in [38] which was 
later updated in [37]. 

References 

[1] BABAR collaboration, B. Aubert et al., Evidence for D^-D" mixing. Phys. Rev. Lett. 98 (2007) 211802; 
BELLE collaboration, M. Staric et al.. Evidence for D° - Z)° mixing^ Phys. Rev. Lett. 98 (2007) 211803; 
BELLE collaboration, L. M. Zhang et. al.. Measurement of — & mixing in D KsTt~^tt^ decays. 
Phys. Rev. Lett. 99 (2007) 131803. 

[2] HFAG Collaboration, D. Asner et. al.. Averages of b-hadron, c-hadron, and r-lepton Properties. 
arXiv:1010.1589. 

[3] L I. Bigi, in: "Heavy Flavour Physics: a Probe of Nature's Grand Design" , Proc. of the International 
School of Physics "Enrico Fermi", Course CXXXVII, lOS Press, (1998) p. 645; arXiv:hep-ph/9712475, 
p. 56/57. 

[4] I. I. Bigi, N. G. Urahsev, D° - L»° Oscillations as a Probe of Quark-Hadron Duality. Nucl. Phys. B 592 
(2001) 92 [hep-ph/0005089]; S. Bianco et al., A Cicerone for the Physics of Charm. Riv. Nuovo Cim., 
26 N7-8 (2003) p. 151 - 153 [arXiv:hep-ex/0309021]. 

[5] M. Bobrowski et. al., How large can the SM contribution to CP violation in — mixing be?. 
JHEP 03 (2010) 009; A. Lenz and M. Bobrowski, Standard Model Predictions for -oscillations and 
CP-violation. arXiv:1011.5608. 

[6] A. Falk et al.,S'f7(3) Breaking and DO-DObar Mixing. Phys. Rev. D 65 (2002) 054034 [arXiv:hep- 
ph/0110317]; The D° - D° mass difference from a dispersion relation. Phys. Rev. D 69 (2004) 114021 
[arXiv:hep-ph/0402204]. 

[7] K. Blum et al.. Combining — Mixing and — Mixing to Constrain the Flavor Structure of 
New Physics. Phys. Rev. Lett. 102 (2009) 211802. 

[8] N. Arkani-Hamed, A. G. Cohen, and H. Georgi, (De) constructing dimensions. Phys. Rev. Lett. 86 

(2001) 47574761, [arXiv:hep-th/0104005]; N. Arkani-Hamed, A. G. Cohen, and H. Georgi, Electroweak 
symmetry breaking from dimensional deconstruction. Phys. Lett. B 513 (2001) 232240, [arXiv:hep- 
ph/0105239]; N. Arkani-Hamed et al.. The minimal moose for a little Higgs. JHEP 08 (2002) 21; N. 
Arkani-Hamed et al., Phenomenology of Electroweak Symmetry Breaking from Theory Space. JHEP 08 

(2002) 020. 

[9] D. E. Kaplan and M. Schmaltz, The little Higgs from a simple group. JHEP 10 (2003) 039; M. Schmahz, 
The simplest little Higgs. JHEP 08 (2004) 056. 



22 



REFERENCES 



[10] M. Perelstein, Little Higgs models and their phenomenology. Prog, in Part, and Nucl. Phys. 58 (2007) 
247. 

N. Arkani-Hamed, A. G. Cohen, E. Katz and A. E. Nelson, The littlest Higgs. JHEP 07 (2002) 034. 

T. Han, H. E. Logan, B McElrath and L-T Wang, Phenomenology of the little Higgs model. Phys. Rev. 
D 67 (2003) 095004. 

H.-C. Cheng and I. Low, TeV symmetry and the little hierarchy problem. JHEP 09 (2003) 051, 
[arXiv:hep-ph/0308199] . 

H. -C. Cheng and I. Low, Little hierarchy, little Higgses, and a little symmetry. JHEP 08 (2004) 061, 
[arXiv:hep-ph/0405243]. 

I. Low, T parity and the littlest Higgs. JHEP 10 (2004) 067. 

J. Hubisz and P. Meade. Phenomenology of the littlest Higgs model with T parity. Phys. Rev. D 71 
(2005) 035016. 

I. I. Bigi et ah, CP Violation in D'^ — Z)" oscillations: general considerations and applications to the 
Littlest Higgs model with T parity. JHEP 07 (2009) 097. 

A. Paul, I. I. Bigi and S. Recksiegel, D'^ — > 77 and — > rates on an unlikely impact of the 

littlest Higgs model with T parity. Phys. Rev. D 82 (2010) 094006. 

SuperB Collaboration: B O'Leary et al., SuperB Progress Reports ~ Physics. arXiv:1008:1541; SuperB 
Collaboration: E. Granges et al., SuperB Progress Report - Detector. arXiv: 1007:4241; SuperB Collab- 
oration: M. E. Biagini et al., SuperB Progress Report - Accelerator. arXiv:1009:6178. 

M. Blanke et al.. Littlest Higgs model with T parity confronting the new data on D'^ — mixing. Phys. 
Lett. B 657 (2007) 81. 

G. Burdman, E. Golowieh, J. Hewett and S. Pakvasa, Radiative weak decays of charm mesonsPhys. 
Rev. D 52 (1995) 6383. 

T. Inami and C. Lim, Effects of Superheavy Quarks and Leptons in Low-Energy Weak Processes Kl — >■ 
Hp., K+ ■K+vv and K° o Prog. Theor. Phys. 65 (1) (1981) 297. 

C. Greub et al.. The c — > wy Contribution to Weak radiative Charm Decay. Phys. Lett. B 382 (1996) 
415. 

G. Burdman et al.. Rare Charm Decays in the Standard Model and Beyond. Phys. Rev. D 66 (2002) 
014009. 

S. Fajfer, S. Prelovsek and P. Singer, Rare charm meson decays D — > Pl^l^ and c — > ul^l^ in the 
standard model and the minimal supersymmetric standard model. Phys. Rev. D 64 (2001) 114009. 

A. J. Buras and M. Miinz, Effective Hamiltonian for B — > Xse'^e~ beyond leading logarithms in the 
naive dimensional regularization and 't Hooft-Veltman schemes. Phys. Rev. D 52 (1995) 186. 

S. Fajfer, P. Singer and J. Zupan, The radiative leptonic decays in the standard model and beyond. Eur. 
Phys. J. C 27 (2003) 201. 

B. Grinstein, M.J. Savage and M.B. Wise, B — >• XsC'^e" in the six-quark model. Nucl. Phys. B 319 
(1989) 271. 

S. Fajfer, S. Prelovsek and P. Singer, Resonant and nonresonant contributions to the weak D — > Vl^l~ 
decays. Phys. Rev. D 58 (1998) 094038. 

A. Ali, T. Mannel and T. Morozumi, Forward-backward asymmetry of dilepton angular distribution in 
the decay b^ sl+l'. Phys. Lett. B 273 (1991) 505. 



23 



REFERENCES 



[31] F. Kriiger and L.M. Sehgal, CP violation in the decay B X^e+e". Phys. Rev. D 55 (1997) 2799. 

[32] G. Buchalla, G. Hillier and G. Isidori, Phenomenology of nonstandard Z couplings in exclusive semilep- 
tonic s transitions. Phys. Rev. D 63 (2000) 014015. 

[33] S. Rai Choudhury, Forward-backward asymmetry in B ^ X^e+e". Phys. Rev. D 56 (1997) 6028. 

[34] J. Hubisz, S. J. Lee and G. Paz, The flavor of a little Higgs with T parity. JHEP 06 (2006) 041. 

[35] M. Blanke et al., Another look at the flavour structure of the Littlest Higgs model with T parity. Physics 
Letters B 646 (2007) 253. 

[36] M. Blanke et al., P article- Antiparticle Mixing.... JHEP 12 (2006) 003,[hep-ph/0605214]; M. Blanke 
et al., Correlations between e'/e and rare K decays in the Littlest Higgs model with T parity. JHEP 
06 (2007) 082; M. Blanke and A. J. Buras, A Guide to Flavour Changing Neutral Currents in the 
Littlest Higgs Model with T parity. Acta Phys. Polon. B 38 (2007) 2923,[arXiv:hep-ph/0703117]; A. J. 
Buras and C. Tarantino, Quark and Lepton Flavour Physics in the Littlest Higgs Model with T parity. 
Proceedings of the Workshop CKM2006 [arXiv:hep-ph/q702202] (2007); M. Blanke et al, The Littlest 
Higgs Model with T parity Facing CP-Violation in B^ — B^ Mixing. [arXiv0805.4393v2]; 

[37] M. Blanke et al., FCNC Processes in the Littlest Higgs Model with T parity: an Update. Acta Phys. 
Polon. B 41 (2010) 657. 

[38] M. Blanke et al.. Rare and CP-violating K and B decays in the Littlest Higgs model with T parity. JHEP 
01 (2007) 066. 

[39] T. Feldmann, C. Proniberger and S. Recksiegel, Characterising New Physics Models by Effective Di- 
mensionality of Parameter Space. arXiv:1009.5283. 

[40] S. Fajfer and S. Prelovsek, Effects of littlest Higgs model in rare D meson decays. Phys. Rev, D 73 
(2006) 054026. 

[41] N. Cabbibo, in Proceedings of the XIII*'' International Conference on High Energy Physics, Berkeley, 
California, 1966 (University of California Press, Berkeley, 1967), 29. 

[42] A. J. Buras, M. Misiak, M. Miinz and S. Pokorski, Theoretical Uncertainties and Phenomenological 
Aspects of B ^ Xs-f Decay. Nucl. Phys B 424 (1994) 374. 

[43] T. Goto et al.. Ultraviolet divergences of flavor changing amplitudes in the littlest Higgs model with T 
parity. Phys. Lett. B 670 (2009) 378. 

[44] F. del Aguila et al., Precise limits from lepton flavour violating processes on the Littlest Higgs model 
with T parity. JHEP 01 (2009) 080. 



24 



